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Within the broad class of multiferroics (compounds showing a coexistence of magnetism and ferro- 
electricity), we focus on the subclass of "improper electronic ferroelectrics" , i.e. correlated materials 
where electronic degrees of freedom (such as spin, charge or orbital) drive ferroelectricity. In par- 
ticular, in spin-induced ferroelectrics, there is not only a coexistence of the two intriguing magnetic 
and dipolar orders; rather, there is such an intimate link that one drives the other, suggesting a 
giant magnetoelectric coupling. Via first-principles approaches based on density functional theory, 
we review the microscopic mechanisms at the basis of multiferroicity in several compounds, ranging 
from transition metal oxides to organic multiferroics (MFs) to organic-inorganic hybrids. 



I. INTRODUCTION 

Materials where cooperative phenomena - such as a 
switchable long-range dipolar or magnetic ordering or a 
structural deformation - spontaneously emerge below a 
critical temperature are termed "ferroic"— . Compounds, 
where more than one kind of ferroic order are estab- 
lished, are consequently denoted as "multiferroics"—"—. In 
this Colloquium paper, we will concentrate on materials 
showing a coexistence between (anti-)ferroelectricity and 
(anti-)ferromagnetism. In those systems, many degrees 
of freedom are simultaneously active with competing en- 
ergy scales, which in turn make the response of multi- 
ferroics to external stimuli (electric and magnetic fields, 
pressure, strain, doping, ...) unusually large. Indeed, 
multiferroics offer a wide playground for colossal cross- 
coupled effects to emerge. What is generally meant by 
"cross-coupling" is a physical response not induced by its 
conjugate field. Examples are magnetoelectricity (change 
in magnetic properties induced by electric fields or in fer- 
roelectric (FE) properties induced by magnetic fields), 
piezoelectricity (change in structural properties induced 
by electric fields or in ferroelectric properties induced 
by structural deformation), magnetostriction (change in 
magnetic properties induced by structural deformation 
or in structural properties induced by magnetic fields) 

Due to their multifunctional nature, multiferroics hold 
great potential for future technological applications (such 
as sensors, memories, actuators, switches)^!^. At the sa- 
me time, their complexity poses serious challenges for 
modelling: it calls for an accurate treatment of correlated 
3d- or 4f-electrons and excited states, as well as for a care- 
ful description of the delicate coupling between electronic 
(spin, charge, orbital) degrees of freedom and structural 
distortions and crystal symmetries. Multiferroics there- 
fore constitute one of the most interesting though chal- 
lenging classes in modern materials science. 

Magnetism is likely one of the earliest discovered phe- 
nomenon in condensed matter, with lodestone proper- 
ties already known around 600 BC; ferroelectricity, on 
the other hand, is a property found rather recently, 
with its earliest observation probably going back to 1921 
when Valasek observed an electric hysteresis in Rochelle 




FIG. 1. Pictorial representation of the different microscopic 
mechanisms that can lead to multiferroicity. 



saltS. The combination between ferroelectricity and mag- 
netism, i.e. multiferroicity, was for long time consid- 
ered to be a very rare phenomenon^-. This common be- 
lief mostly derived from an early observatior^ii: most of 
the ferroelectric oxides {i.e. BaTiOs, PbTiOa, LiNbOa) 
have a perovskite-like structure where the perovskite B- 
site cation, which is the one that mostly off-centers and 
gives rise to ferroelectricity, shows an electronic dP con- 
figuration {i.e. Ti^"*", Nb^"*"). The empty d-shell - re- 
quired for ferroelectricity - therefore seemed to preclude 
any coexistence with magnetism. More generally and 
very recently, the the conditions of multiferroicity in d" 
perovskites are derived from the pseudo Jahn- Teller ef- 
fect, due to which ferroelectric displacements are trig- 
gered by vibronic coupling between ground and excited 
electronic states of opposite parity but same spin mul- 
tiplicity; it takes place for some specific d" configura- 
tions and spin states only. In combination with the 
high-spinlow-spin crossover effect this leads to a novel 
phenomenon, the magnetic-ferroelectric (multiferroics) 
crossover which predicts magnetoelectric effects with ex- 
citing functionalities including electric magnetization and 
demagnetizationiiS, 

"Conventional" multiferroics (whose class prototype is 
BiFeOgi^"— ), are materials where, e.g., lone-pair elec- 
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trons of the A-site cation in the perovskite structure 
(such as Bi s with high polarizabihty) gives origin to 
the ferroelectric order, whereas the magnetic order is de- 
termined by the exchange interaction among uncompen- 
sated spins of the B-site cation. When this conventional 
approach is pursued, it follows that magnetism and fer- 
roelectricity come from two different atomic sublattices 
and have a chemically different origin, thereby result- 
ing in different ordering temperatures and presumably 
a small magnetoelectric (ME) coupling (although this is 
not always the casei^). These drawbacks could be over- 
come in materials where magnetism and ferroelectricity 
would share the same microscopic origin (i.e. the two 
phenomena physically originating from the same chemi- 
cal species) and where the expected coupling would there- 
fore be much stronger than in conventional multiferroics. 
Indeed, materials like those attracted an incredible at- 
tention in the last decade and will be the main topic of 
the present review. 

Being multiferroics a rather new field of materials sci- 
ence, in most cases the nature of the coupling between 
ferroelectricity and magnetism is unknown and the mi- 
croscopic mechanisms at its basis have to be discovered. 
Very often, the magnetoelectric coupling is mediated by 
different degrees of freedom and interactions (orbital or- 
dering, spin-orbit coupling, charge disproportionation, 
etc), therefore becoming rather complex to unveil. In 
this context, density functional theoryi^ can be of para- 
mount importance, since this "ab-initio" approach is in 
principle able to describe all the many active degrees of 
freedom within the same level of accuracy, at variance 
with model-Hamiltonian approaches where a choice has 
to be made from the very beginning about the relevant 
interactions to be taken into account. 

The challenging issue of multiferroicity motivated our 
recent theoretical activity, which has been mainly focused 
on investigating microscopic mechanisms that could lead 
to ferroelectricity induced by breaking of inversion sym- 
metry {i.e. the necessary symmetry condition for ferro- 
electricity to develop) via electronic degrees of freedom, 
such as charge, spin, orbital orderi^i^. In these latter 
cases, in the framework of phase transitions, ferroelec- 
tricity is said to be ^Hmproper'\ meaning that ferroelec- 
tric polarization appears as a "secondary" order param- 
eter driven by a "primary" electronic order parameter. 
An example will make things clearer: if one considers a 
(primary) antiferromagnetic (AFM) phase transition in 
which the AFM ordering pattern breaks inversion sym- 
metry, then polarization (secondary order parameter) is 
allowed only when the AFM order (primary order pa- 
rameter) sets in. This is what happens, for example, in 
ortho-HoMnO^il^, TbMnsOgiS, CasCoMnOf^. Inci- 
dentally, we here remark that conventional ferroelectric 
oxides show "proper" ferroelectricity: polarization is the 
only order parameter in the phase transition and is not 
driven by (nor drives) any other electronic phase transi- 
tion. 

In what will follow, we will generally consider improper 



ferroelectricity as electronic ferroelectricity, i.e. driven 
by a spin, charge or orbital phase transition. In this 
regard, we note that improper ferroelectricity can oc- 
cur even when the primary order parameter is a (non- 
polar) structural distortion (as already well investigated 
in Refsi^i), and not necessarily an electronic order pa- 
rameter. 

Further related to electronic ferroelectricity, we'd like 
to make the following comment. Polarization can be de- 
fined as the sum of an electronic and of an ionic contribu- 
tion, which can be roughly understood as resulting from 
the polar charge rearrangement and polar distortions, re- 
spectively. Indeed, the definition of polarization is rather 
tricky, in particular in a periodic system. We'll here only 
recall some of the relevant points and we refer the inter- 
ested reader to Refsi^"— where the topic is treated in a 
deeply detailed way. A rough way of defining polarization 
is through the so called "point-charge model" (PCM), 
obtained as the sum of the displacements of the atoms 
(with respect to a reference centrosymmetric paraelectric 
(PE) structure), each multiplied by the nominal valence 
of the corresponding ion. In this case, one adopts a purely 
ionic model and therefore neglects the details of the real 
charge distribution resulting from hybridization and co- 
valent effects. A rigorous way to overcome the difficulties 
in defining polarization is therefore provided in Refsi^"— 
by two equivalent quantum-mechanical approaches, one 
based on the Berry phases and one based on Wannier 
function centers. According to the latter, the electronic 
charge is considered to be localized at the Wannier func- 
tion centers, while the ionic charges reside on the nuclear 
positions, so that the change in electronic polarization 
can be obtained as a vectorial sum of the displacement of 
the Wannier-function centers. Following this approach, 
it is in principle possible to have an electronic contri- 
bution to polarization, due exclusively to Wannier cen- 
ter displacements (say, following a non-centrosymmetric 
AFM arrangement) on top of a centrosymmetric atomic 
configuration (i.e. that would otherwise lead to a vanish- 
ing contribution within simpler classical models, such as 
PCM). We however remark that there cannot be, strictly 
speaking, a purely electronic ferroelectricity, since polar 
atomic displacements induced by a polar electronic order 
are always present, no matter how small. Furthermore, 
one of the information that can actually be provided by 
a first-principles approach used by the present authors 
is indeed the quantification of electronic vs ionic contri- 
bution to polarization, as discussed below. In summary, 
what we mean in this Colloquium by "electronic ferro- 
electricity" is the phenomenon by which a primary elec- 
tronic phase transition breaks inversion symmetry and 
gives, as a by-product, a ferroelectric polarization. 

A general characteristic of electronic ferroelectricity is 
that the expected magnitude of the polarization is likely 
smaller than in standard ferroelectrics; in the case of 
"improper" ferroelectricity, in fact, the polarization is 
driven by a polar electronic charge rearrangement, ac- 
cording to which - as a by-product - the atoms are slightly 
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displaced in a non-centrosymmetric way. On the other 
hand, in the latter "proper" case {i.e. BiFeOa), the in- 
version symmetry is broken primarily by structural dis- 
tortions, involving ionic displacements of the order of a 
tenth of an Angstrom (10 or 100 times larger than in elec- 
tronic improper ferroelectrics). However, electronic fer- 
roelectrics might offer a significant advantage over stan- 
dard ferroelectrics, as for what concerns the switching 
time-scale and the so called "fatigue". In general, the 
latter term means the deterioration of the hysteresis loop 
and the decrease of switching charge after many polariza- 
tion reversal cycles. During switching processes in proper 
ferroelectrics, it is the (slow and heavy) ions that dis- 
place, whereas in improper electronic ferroelectrics it is 
the (fast and light) electrons that move: this fundamen- 
tal difference is definitely expected to lead to a much 
quicker switching time^"— and to a strong reduction of 
fatigue-related problems (although the origin of the latter 
is presently not well understood). 

As for the relevant materials, the natural class where 
one expects multiferroicity to arise is represented by tran- 
sition metal oxides, since they are probably the richest 
class in materials science where structural and electronic 
degrees of freedom are all simultaneously active and in- 
teracting. However, over the years, we showed that nice 
effects can be found in organic materials as wel P^'^^ , 
where various correlation phenomena usually show up 
and, additionally, low dimensionality can play a relevant 
role and can be considered as an additional degree of 
freedom to be tuned and exploited. Finally, very re- 
cently we discovered multiferroicity in organic-inorganic 
hybridsSS, such as perovskite-based metal-organic frame- 
works, where the unlimited variety of organic functional 
groups is nicely joined to the rich functionality arising 
from the perovskite network. 

In this Colloquium paper, we put the emphasis on the 
rich collection of microscopic mechanisms that can lead 
to "electronic ferroelectricity" , each of them manifestly 
at play in a different material and highlighted by means 
of state-of-the-art first-principles calculations in the dif- 
ferent Sections. A pictorial representation of the differ- 
ent possible origins of multiferrocity is shown in Figure 
[TJ In particular, following a Section with computational 
technicalities (see Section [H]), we will first discuss how 
the spin ordering can break inversion symmetry (see Sec- 
tion |lirB]), by introducing general mechanisms based on 
the Heisenberg symmetric exchange and Dzyaloshinskii- 
Moriya antisymmetric exchange. Specific examples will 
be discussed, based on the Heisenberg exchange involv- 
ing exchange coupling between 4/ and 3d electrons in 
orthoferrites (see Section UlI B ip and V-V exchange cou- 
pling in vanadium-based spinel (see Section fill B 2[) . In 
Sec lIII CI we will discuss how a specific charge-ordering 
can give rise to a polarization, either combined with spin- 
ordering or by itself, as extensively reported for mag- 
netite in Sec lIII CTl Orbital order can also in princi- 
ple break inversion symmetry and induce ferroelectricity; 
however, we are not presently aware of any material in 



which 00 is the primary and unique cause of polarization 
(with the possible exception of double-layer manganite, 
although the origin of ferroelectricity is still not well un- 
derstood in that compound). We argue, anyway, that 
00 can be an important ingredient which cooperates 
with hydrogen bonding network to induce ferroelectricity 
in a Cu-based metal-organic framework, as reported in 
Sec lIII D"T1 The case of a donor-acceptor molecular crys- 
tal, such as TTF-CA (see Section IlII El) is an example of 
cooperation between charge transfer, structural dimeriza- 
tion and possibly spin-Peierls transition, that ultimately 
leads to a paradigmatic organic ferroelectric crystal. In 
Sec lIII Fl we focus on recent developments in the field, as 
we'll discuss the case of several multiple non-polar insta- 
bilities which finally result in ferroelectricity, as shown 
by our ab-initio calculations for NaLaMnWOg double- 
perovskites. In Sec lIVI we draw some conclusions and 
offer a perspective view of the field in the near future. 



II. GENERAL COMPUTATIONAL 
FRAMEWORK 

A. Density functional: technicalities 

All the calculations presented here have been per- 
formed within Density Functional Theory (DFT) using 
the Vienna- ab-initio Simulation Package (VASP )^^'^^ . 
The electron-ion interaction is described by PAW 
potentials^^i^ using a plane-wave basis set with appro- 
priate energy cut-off. For Brillouin zone integrations we 
used the Monkhorst-Pack schemes. We refer to the orig- 
inal papers (each mentioned in the appropriate section 
below) for the specific values and details of the compu- 
tational parameters for the different compounds. 

For the exchange-correlation functional, E^^c, we used 
several approximations depending on the specific prob- 
lem at hand. It is well known that there are wide classes 
of materials where density-functional methods fail not 
only quantitatively but also qualitatively. Typical situa- 
tions are materials with localized orbitals, e.g. of tran- 
sition metals or rare earths. A second important class 
is comprised of the correlated organic crystals like low- 
dimensional organic charge-transfer salts. 

The localized nature of the 3d electronic states limits 
to the applicability of common density-functional meth- 
ods like the local density approximation (LDA) or gen- 
eralized gradient approximation (GGA). In fact, these 
standard approximations introduce a spurious Coulomb 
interaction of the electron with its own charge, i.e., the 
electrostatic self-interaction is not entirely compensated. 
This causes fairly large errors for localized states (e.g., 
Mn d states). It tends to destabilize the orbitals and 
decreases their binding energy, leading to an overde- 
localization of the charge density^. The most com- 
monly applied Generalized-Gradient Approximation, the 
Perdew-Burke-Ernzerhof parametrization (PBE), is of- 
ten insufficient to treat with these problems. To over- 
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come these failures of standard DFT approaches, we 
make use of possible ways out. One common approach 
is the DFT+U method^^"-^, where a Hubbard-like U 
term is introduced into the DFT energy functional in 
order to take correlations partially into account. The 
method usually improves the electronic-structure de- 
scription, but it suffers from shortcomings associated 
with the U-dependence of the calculated properties^^Ti^. 
Unfortunately, there is usually no obvious choice of the U 
value to be adopted; common choices are usually based 
either on experimental input or are derived from con- 
strained DFT calculations^. 

Another approach which is becoming widely used in 
the solid-state community is the use of hybrid func- 
tionals. Hybrid functionals go beyond the usual Kohn- 
Sham formalism and fall within the generalized Kohn- 
Sham realization of DFT-^. These functionals consider a 
weighted mixture between exchange defined in Hartree- 
Fock theory using DFT Kohn-Sham orbitals with DFT 
exchange. Popular hybrid functionals B3LYP— , PBEO^, 
and HSE'*^ have been constructed to give good structural, 
thermodynamic, and bonding properties of solids^i^. In 
particular, we used the HSE06 functional, which is par- 
ticularly suitable for solid state applications''^"'*^. 

The ferroelectric polarization is calculated within 
the modern theory of polarization22i2^. where one 
computes the difference of electric polarization, i.e. 
AP^Pfe~Ppe 

=APjo„-|- APgfe, where the subscripts FE, PE, ion and 
ele denote ferroelectric, paraelectric, ionic and electronic 
contribution, respectively. Pion is calculated by summing 
the position of each ion in the unit cell times the number 
of its valence electrons. The electronic contribution is 
obtained by using the Berry phase formalism^'^'^'* . 



B. Symmetry analysis 

The use of symmetry analysis in this field is very rel- 
evant, as will be evident in the following sections. We 
recall in this framework that, from the symmetry point 
of view, a magnetic ordering breaks time-reversal sym- 
metry, a ferroelectric ordering breaks inversion symme- 
try, so that both time and space inversion symmetry are 
absent in multiferroics (although by two different order 
parameters, such as for example magnetization and po- 
larization). Symmetry analysis can be very helpful in 
identifying polar ionic or electronic arrangements or in 
suggesting a physically-sound paraelectric structure {i.e. 
reasonably close to the ferroelectric structure, in terms 
of ionic displacements), that can be taken as reference 
paraelectric structure when calculating the polarization. 
Furthermore, useful tools for the calculation of the po- 
larization as well as for its analysis, are those of the Bil- 
bao Crystallographic Serveri^rSi^ or of the ISODISTORT 
web site^. 



III. MECHANISMS AND MATERIALS 

A. Introduction: magnetic interactions and spin 
Hamiltonian approach 

In general, the interactions between magnetic centers 
in condensed matter systems as well as in molecular 
systems, have a twofold orgin, one is purely magnetic 
and the other elettrostatic in nature and they can be 
described through-space and through-bond, respectively. 
The former is the usual point dipolar approximation be- 
tween magnetic centers, which can be safely disregarded 
in most cases. The latter relies on the electrostatic in- 
teraction responsible of the formation of the chemical 
bonds. The states of the the interacting centres are de- 
scribed by a set of orbitals. As a rule of thumb for the 
magnetic interaction, if the single occupied orbitals are 
orthogonal to each other, the two spins of the electrons 
will be parallel to each other (ferromagnetic coupling), 
whereas if the orbitals have a non-zero overlap the spins 
will tend to orient antiparallel to each other (antiferro- 
magnetic coupling)^^'^^. 

From an historical point of view, the description of 
the magnetic interaction was performed by using local- 
ized magnetic orbitals or a valence-bond approach^"—. 
Other approaches have been developed since then, which 
rely on the tight-binding approaches j^i*^ or density func- 
tional theory. ^■^"^^ The basis of the magnetic interaction 
is the antisymmetric nature of the total wavefunction, 
which shows up as an effective exchange interaction. The 
exchange interaction may occur directly (direct interac- 
tion) or through a formally diamagnetic ligand (super- 
exchange). The famous Goodenough-Kanamori rules 
represent a qualitatively account of the features of the 
magnetic coupling between different centers. 

It is often useful to introduce a spin Hamiltonian in 
order to eliminate all the orbital degree of freedoms and 
replace them with spin coordinates. A central approxi- 
mation for such a mapping is that orbital moment is es- 
sentially quenched as it often occurs in solids, and it can 
be eventually treated as perturbation. The spin Hamil- 
tonian approach can be used for treating: 

• Zeeman and crystal field terms for isolated ions; 

• electron nucleus interaction terms (hyperfine inter- 
actions); 

• interaction betweeen spin pairs. 

Hereafter, we will focus on the last term. It can 
be represented be a spin-spin Hamiltonian, which 
can be written as H=S'i • J12 • S2 where Si^2 are 
the spin operators for the magnetic center 1 and 2, 
respectively. J 12 is a matrix which describes the in- 
teraction, which in general is not symmetric and may 
have a non-zero trace. It is always possible to break- 
down the J12 tensor equivalenty into three contributions: 
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-Jl25i • 52 + Sv Di2 ■S2+di2 -01 X 52) 

where Ji2=-(l/3)Tr J12; D^/ = (l/2)(J5^f )+ jf^"- 

S^p (l/3)Tr(Ji2); di2=(l/2) (jf^^ - J^f ) and a,/3,7 are 
Cartesian components. The first term in the previous for- 
mula is referred to as the isotropic which tend to keep the 
spins coUinear; the second as the anisotropic and it tends 
to orient the spins along a given orientation in space; the 
third as the antisymmetric spin-spin contribution to the 
magnetic interaction and it tends to cant them by 90°. 
In many cases, the first term is the dominant one, and 
the other terms can be introduced as perturbation. Here, 
for positive J12 we have ferromagnetic coupling, while for 
antiferromagnetic coupling J12 is negative. The mecha- 
nisms at the basis of the exchange interactiona have been 
first introduced and discussed by Anderson, and traslated 
in rulc-of-thumbs by Goodenough and Kanamori. 

In the following, we want to discuss the origin of the 
phenomenon of weak ferromagnetism and how it arises 
from the previous spin Hamiltonian. The origin of weak- 
ferromagnetism is usaully ascribed to the so-called 
Dzyaloshinskii-Moriya (DM) interaction. This exchange 
interaction results from the interplay of the Coulomb in- 
teraction and the spin-orbit coupling in systems of low 
crystal symmetry. The DM interaction is an important 
term playing a crucial role in may magnetic systems. The 
weak ferromagnetism is characterized by a small net mag- 
netic moment resulting from spin moments that nearly 
cancel each other. It was first observed in haematite, 
a— Fe203''^'^^.Dzialoshinki showed that it was an in- 
trinsic effect due to the particular symmetry proper- 
ties of the crystal structure and the magnetic moments 
arrangments2S. Furthermore, it showed that second term 
and the third term in the previous spin Hamiltonian, 
i.e. magnetocrystalline anisotropy (or anisotropic term) 
and the anisotropic exchange (or antisymmetric spin- 
spin interaction) rispectively, can lead to a small fer- 
romagnetic moment in an otherwise antiferromagnetic 
crystal. Moriya^° showed that the Dzialoshinski's ex- 
planation can be interpreted in the framework of Ander- 
son's perturbation approach to magnetic superexchange. 
Furthermore, he showed that depending on the type of 
crystal structure either of two mechanism, magnetocrys- 
talling anisotropy or antisymmetric exchange, can be the 
source for the canting of magnetic moments. For ex- 
ample, for a— Fe203, it is the antisymmetric exchange 
that plays the dominant role whereas, in the case of 
NiF2 , antisymmetric exchange is ruled out in favor of the 
magnetocrystalling anisotropy which is here the impor- 
tant term giving rise to the weak ferromagnetic compo- 
nent. In the triangular antiferromagnetc MnaSn the an- 
tisymmetric exchange contributions from different atoms 
cancel perfectly and can not be the reason for the ob- 
served weak ferromagnetism. Here, the magnetocrys- 
talline anisotropy term is the source of the spin cant- 
ing. Another example, is the Cu based metal-organic 
frameworks, where it can be shown that antisymmetric 
exchange is zero for symmetry, while the only term con- 



tributing to the spin canting is the magnetocrystalline 
anisotropy-Si. 

B. Spin ordering 



a) AFM-E, +P 



b) AFM-E, -P 



c) CW, +P 



d) CCW, -P 




FIG. 2. Examples of polar spin configurations (shown in the 
Mn02 basal plane of Prima manganites) : a) and b) coUinear 
cases (such as AFM-E in HoMnOa) with opposite polariza- 
tion, induced by flipping the direction of two spins, c) and 
d) spiral cases (such as TbMnOa) with opposite polarization, 
induced by the different clockwise (CW) vs counterclockwise 
(CCW) rotation of spins in the vertical plane. The conse- 
quent exchange-strictive mechanisms are also shown: semi- 
transparent oxygens schematically show the displaced atoms 
(exaggerated for clarity) occurring upon spin ordering: In 
a) and d) the spin configuration is such that oxygens move 
"down", whereas in b) and c) oxygens move "up". 

In recent years, people have focused mainly on two dif- 
ferent mechanisms of magnetically-induced ferroelectric 
polarization, P, driven by either symmetric Si • Sj Heisen- 
berg (H) exchang o^^'^^ or antisymmetric Si x Sj spin- 
orbit-like Dzyaloshinskii-Moriya (DM) exchange^"—. 
The prototypical materials taken as representative ex- 
amples of the two mechanisms are by far HoMnOa and 
TbMnOa, respectively: i) the first manganite shows 
a coUinear t ~ t ^ i ^ i spin configuration along 
[101] and [10-1] directions in the Pnma setting ; al- 
ternatively, this peculiar AFM ordering can be seen as 
zig-zag ferromagnetic chains antiferromagnetically cou- 
pled in the ac Mn02 plane . This so-called AFM-E 
spin order (which is stable when the A-site ionic radii 
is small, such as Tm, Lu, Yb rare-earth cations in ortho- 
manganites) clearly lacks inversion symmetry (see Figure 
[5] a). From the microscopic point of view, the polar- 
ization comes from the inequivalency between the oxy- 
gens bonded to two Mn with parallel spins, with respect 
to those linked to Mn with antiparallel spins, as dis- 
cussed in detail in Refsj^iiii^ ii) the second manganite, 
TbMnOa shows as ground-state a spin cycloidal spiral 
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in the be plane. In that case, non-cohinear spins induce 
polarization, as predicted by the formula P (x en x Q, 
where Q is the spiral wave- vector parallel to the chain 
direction and e„ oc S'„ x 5'„+i is the spin-rotation axis 
[see Figure [2] c) and d)] The vector product Sn x Sn+i 
can be shown to be proportional to the spin current jg, 
in turn linked, via a vector-potential relation^^, to the 
Dzyaloshinskii-Moriya interaction, so that the mecha- 
nism of ferroelectricity in spin-spirals is generally labeled 
as "spin-current-induced"—. 

In the context of spin-induced ferroelectricity, it is in- 
teresting to discuss how "switching" of polarization oc- 
curs. We recall that, in standard ferroelectrics, ionic dis- 
placements occur in equal magnitudes and opposite di- 
rections (with reference to the paraelectric state) when 
reaching the "-|-P" or "-P" state. However, in spin- 
induced ferroelectricity, since polarization is determined 
primarily by the spin ordering and not by ionic displace- 
ments, one expects some changes to occur in the ar- 
rangement of magnetic moments. For example, in AFM- 
E ortho-manganites, switching occurs when half of the 
spins in the unit cells have their direction flipped by 180°, 
so that all oxygens that were previously connected to par- 
allel spins now become oxygens being bonded to antipar- 
allel spins (and viceversa, cfr Figure [5] a) and b)). When 
focusing on spiral TbMnOa, switching occurs when re- 
versing the vector spin chirality, i.e. when reversing a 
clockwise spiral into an anticlockwise spiral, as shown in 
Figure [2] c) and d). 

What is important to emphasize for spin-driven ferro- 
electricity is again the difference between the electronic 
and ionic contributions to polarization: the first arises 
from the polar rearrangement of the electronic charge, it 
is present even with ions arranged in a centrosymmet- 
ric way and was shown to be relevant in both the spi- 
ral and AFM-E cases; the second is mainly referred to 
as exchange-striction (symmetric for the Heisenberg case 
and antisymmetric for the DM case), meaning that the 
ions move, when the spin ordering is established, to gain 
energy from exchange terms (see Figure[2]for a schematic 
representation) . Which of the two - electronic vs ionic - 
contributions is more relevant actually depends pretty 
much on the considered system: from first-principles cal- 
culations, the two were found to be almost the same in 
AFM-E HoMnOjii, whereas the ionic was estimated to 
be much larger in spiral TbMnOg^. 

When comparing the two DM and H mechanisms, one 
expects that, being the former DM-induced polarization 
driven by relativistic effects which are not so strong in 
3d-based materials, the DM-related polarization should 
be weaker than the H-induced one. Indeed, in a variety 
of materials (ranging from nickelates, such as LuNiOs*^, 
to manganites, such as HoMnOg^^, or sulfides, such as 
Cu2MnSnS4^), we have shown that the size of P in- 
duced by the relativistic DM interaction (occurring in 
spin-spiral-based oxides) is much lower than that caused 
by the H interaction. Related to this, a comment is in or- 
der: it was very controversial how large the polarization 



actually was in AFM-E like rare-earth ortho-manganites, 
a paradigmatic case of Heisenberg-driven ferroelectricity. 
Earlier experiments from Lorenz et alS^ reported ferro- 
electricity on polycrystalline samples of ortho- HoMnOa 
of the order of 10^"^ /iC/cm^, three orders of magnitude 
smaller than what predicted by our DFT calculations. 
Recently, however, the group of Tokura,^*' by means of 
advanced growth techniques and measurements, focused 
on several rare-earth ortho-manganites (RMnOa, R — 
Ho, Lu, EuxYi-x) estimated the genuine values of P 
(about 0.5 /xC/m^) in the E-type phase, which is more 
than 10 times as large as that of the be cycloidal phase. 
Slightly earlier, a theoretical model was reported in the 
context of electromagnon excitations in RMnOa^^. One 
of the outcome was the estimate of the polarization in E- 
type manganites based on optical absorption data mea- 
sured for TbMnOa in the spiral-phase: P was found to 
be of the order of 1 /xC/cm^. Moreover, Pomyakushin 
et al^ have reported the polarization of about 0.15 
/iC/cm^ for E-type TmMnOa. There is therefore now 
a growing consensus on the possibility of breaking inver- 
sion symmetry (therefore paving the way to improper 
ferroelectricity) via the symmetric magnetic exchange 
(Heisenberg-like) in coUinear frustrated systems, in addi- 
tion to the (well consolidated) analogous effect in the an- 
tisymmetric counterpart (Dzyaloshinskii-Moriya-like) of 
magnetic exchange. 

In what follows, we'll discuss two peculiar cases in 
which the Heisenberg exchange striction is at play to in- 
duce ferroelectricity, according to a microscopic mecha- 
nism based on the up- up- down- down spin chain: a rare- 
earth orthoferrite and Cd-based vanadate. Despite the 
similarity with the prototypical HoMnOa, both systems 
show some peculiarities: in the first case, it is the inter- 
action between 4/ and Zd spins, a-priori not expected to 
be very strong, that causes a sizeable P, whereas in the 
second case, it is the peculiar spinel coordination - and 
related oxygen arrangement - which makes the mecha- 
nism underlying ferroelectricity more intriguing. 



1. f — d coupling in DyFeOs 

DyFeOa has been studied since a long time for 
its interesting physical properties^^^-i^^. Dysprosium- 
Orthoferrite, DyFeOa (cfr Figure [3] a), was suggested 
from 

experiments^*^^ to show, below the ordering temperature 
of Dy and upon application of a magnetic field parallel 
to the c-axis and larger than a specified critical mag- 
netic field, a MF phase, with weak ferromagnetism and 
ferroelectricity (M'^0.5 /is per formula-unit and P ^ 0.2 
/iC/cm^). By performing accurate DFT calculations (see 
details in Refi^Sl), we focused on the mechanism at the 
basis of the polar behavior, neglecting the weak magnetic 
moment that most likely arises from a Fe spin-canting 
(not directly related to ferroelectricity) and considering 
a simpler collinear spin arrangement. We showed that 
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a) Pnma structure b) Spin configuration 




FIG. 3. Orthorhombic DyFeOs: a) crystal structure (FeOs 
octahedra and atomic planes highlighted); b) spin config- 
uration: blue (red) arrows denote Dy (Fe) spins. Black 
bold arrows denote the displacement (with respect to the 
non-magnetic paraelectric structure) of Fe02 planes, due to 
exchange-striction. 



4/ electrons, more-often-than-not neglected in modelling 
studies, play an unexpected and important role in sta- 
bilizing the magnetic-field-induced ferroelectric state of 
DyFeOs. Indeed, we found the FE polarization to be 
mainly driven by an exchange-strictive mechanism, work- 
ing between adjacent spin-polarized Fe and Dy layers ar- 
ranged in a Fet-Dyt-Fe4,-Dy4, fashion (see Figure [3] b) : 
DyO atomic planes move towards Fe02 ionic planes so as 
to maximize the energy gain coming from a Fe-3d/Dy-4/ 
ferromagnetic coupling. Here, at variance with AFM-E 
HoMnOa where the ttii spins are all equivalent, we have 
an intrinsic inequi valency between Fe-d and Dy-/ spins, 
so that Fet-Dyt-Fe4,-Dy4, chains "dimerize" and give rise 
to polarization. Indeed, the distortions lead to an al- 
ternate short-long-short-long interlayer distance between 
DyO and Fe02 planes, with dpM = 1.898 A and dAFM 
= 1.910 A, to be compared with dideai = 1.904 A in the 
unrelaxed non-spin-polarized case (cfr Figure[3]). Indeed, 
we were able to identify two degenerate and switchable 
polar states, i.e. characterized by a sign-reversal of the 
FE polarization (±P) and connected by a relative ro- 
tation of the direction of Dy spins (with respect to Fe 
spins). The estimated FE polarization, in good agree- 
ment with experimentsiSS. shows an unexpectedly large 
magnitude (~0. 1-0.2 ^C/cm^). 

So far we have discussed the exchange striction mech- 
anism in the Heisenberg framework. Obviously, this cou- 
pling must have a microscopic explanation in terms of 
interaction between orbitals. In order to unveil the mi- 
croscopic origin of the polarization, a careful analysis of 
atomic displacements as well as of electronic structure is 
needed. 

Let's consider the virtual paraelectric phase where 



the ions are locked at a centrosymmetric (OS) position 
{Pnma setting, point group D2h)- A suitable spin con- 
figuration which does not break the inversion symmetry 
shows the Dy spin intralayer FM coupled, but rotated 
with respect to the Fe spins by 90°. The FEi (+P) or 
FE2 (— -P) state can be obtained by progressively rotat- 
ing the Dy spins in the ac plane, counter-clockwise for 
FEi or clockwise for FE2, when viewed from the posi- 
tive direction of the b axiaiSl. In this way, one recovers 
the AFM-A spin configuration, i.e. ferromagnetic planes 
antiferromagnetically coupled along the b axis. 

From symmetry point of view, the configuration with 
ortogonal spins has the magnetic space group P212121 
(N. 19.1.119), and therefore the space group for the nu- 
clear sites is P212121 (N. 19). This means that the mag- 
netic ordering breaks all the symmetry operations con- 
taining time inversion, plus the inversion centre and all 
mirror planes. However, we remark that, although the 
magnetic ordering breaks the inversion centre, the re- 
sulting space group symmetry is non-polar. It should be 
noted that the nuclei are not constrained by this symme- 
try to stay in the ideal Pnma configuration that we used, 
and can in principle relax to a more general arrangement 
compatible with the P212121 space group, by means of a 
non-polar distortion. Also the magnetic symmetry does 
not force the magnetic moments to be strictly orthogo- 
nal. Both Fe and Dy magnetic moments could have some 
antiferromagnetic components along some of the other 
axes (with different sign correlations among sites In 
our work, we have disregarded these (presumably small) 
deviations from collinear spin configuration. 

The configuration with collinear spins has the mag- 
netic space group Pn2\a (N. 33.1.226 in non-standard 
setting), and therefore the space group for the nuclear 
sites is Pn21a (N. 33 in non standard setting). This 
means that the magnetic ordering breaks all the symme- 
try operations containing time inversion, plus the inver- 
sion centre, the binary axes along x and z, and the mirror 
plane perpendicular to y. This magnetic ordering breaks 
the inversion centre as in the previous case, but now the 
resulting space group symmetry is indeed polar along the 
6-axis. Again, the nuclei are not constrained to stay in 
the ideal Pnma configuration that we used, and could re- 
lax to a more general arrangement compatible with the 
space group Pn21a, through a distortion which will be 
polar along y. So in principle, it is quite similar to the 
first case above, the difference being that the possible dis- 
tortions relaxing the nuclear positions will be polar, and 
therefore can yield some macroscopic polarization along 
the 6-axis, while in the previous case the possible relax- 
ations of the nuclear structure are necessarily non-polar. 
As in the previous case also, the magnetic symmetry of 
the configuration does not force the magnetic moments 
to be strictly collinear, and some additional antiferro- 
magnetic arrangements of other additional components 
of the magnetic moments (with different sign relations 
among the sites) are allowed by symmetry, and, however 
small, they will in principle be present in a fully relaxed 
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structure^"^. Also in this case, we have not considered 
these spin components. 

Let's focus on equatorial oxygens, Oeq, which occupy 
the 8d Wyckoff positions (WPs). In the FE phase, when 
the symmetry is lowered to C2V, the Of,q become inequiv- 
alent and a WP splitting 8d —>■ Aa + 4a shows up. In- 
spection into the local spin configuration around OeqS 
explains the reason of this inequi valency: the O^q sand- 
wiched by FM coupled Fe and Dy layers, have two tFe 
and two tDy atoms as nearest neighbors (we call them 
as Ol'^)] when sandwiched by Fe and Dy layers AFM 
coupled, they have two tFe and two |Dy atoms as near- 
est neighbors (labelled as 0|^^). The inequivalency due 
to the local spin environment is confirmed by the fol- 
lowing computational experiment: if we impose the FEi 
(or FE2) spin configuration on top of the centrosymmet- 
ric (CS) ionic structure, 0|^^ and 0|g^ become inequiv- 
alent: 0|^^ has ±0.194 fiB and 0|^^ has ±0.207 as 
induced spin moment. This time, however, no WP split- 
ting is involved, since the ions are frozen in CS positions. 
To rule out any numerical artifact on this small differ- 
ence, we impose the PE spin configuration on top of the 
CS ionic structure. In this case, all OeqS carry induced 
spin moments of exactly the same magnitude, becoming 
equivalent again. This leads to the conclusion that the 
change of spin state in going from PE to FEi must be 
the source of the inequivalency of OeqS, and, eventually, 
it should be strongly correlated to the presence of ferro- 
electricity. If so, the ferroelectric state in our toy-model 
is spin-induced. To support this conclusion we note that 
in the PE spin configuration on top of the CS positions, 
Ptot=0 while in FEi, Ptot is different from zero, even 
when the ions are at CS positions. 

In passing we note that all oxygens remain equivalent 
when the Dy-/ electrons are not treated (as done so far) 
as valence electrons, but they are treated as "frozen" 
in the core. In this computational experiment, the Dy 
atoms lose their spins and Ol'^—>- Ot.nospm g^^^ 0|^-^-> 
Qi,«ospm. Qt,nosp»« equivalent to 0^'"°*^™ since here 
we are neglecting the spin-orbit coupling. We are, there- 
fore, led to the conclusion that a signature of the FE in- 
stability is the spin- induced inequivalency of Oeq, which, 
in turn, must be correlated to Dy-/ states, which carry 
the Dy spins. An analysis of the symmetry breaking 
distortionSi^^— sheds further light into the microscopic 
mechanism. The mode decomposition confirms that a 
polar mode is involved, called GM4— . The correspond- 
ing pattern of atomic displacements (not shown here, for 
details see RefjiSZ) with respect to the CS structure high- 
lights the subtle inequivalency of OeqS: Ol'^ (0|^^) move 
in such a way to decrease (increase) the distance to its 
neighbor Dy atom. For 0|^^, doy-o is 2.478 A; for Ol--^, 
doy-o is 2.496 A (the corresponding distance in the PE 
phase is 2.487 A). This would suggest that a weak bond- 
ing interaction is active between the FM layers, in turn 
responsible for the changes in the distances the magnetic 
layers, i.e. for the dimerization and the rising of the po- 



larization. 

A useful tool to study tiny differences in bonding inter- 
action in solid state systems is the electron localization 
function (ELF)iiSiiii. The electron localization function 
was introduced by Becke and Edgecombe as a measure 
of the probability of finding an electron in the neighbor- 
hood of another electron with the same spinJ^iiii ELF 
is thus a measure of the Pauli repulsion, which is active 
when same spin-electron wavefunctions start to overlap 
causing a repulsion due to antisymmetrization postulate. 
The ELF values lie by definition between zero and one. 
Values are close to 1, if in the vicinity of one electron 
no other electron with the same spin may be found, for 
instance as occurring in bonding pairs or lone pairsW. 

We here look for a signature that two ferromagnetic 
layers are interacting through exchange striction. In 
terms of chemical bond picture, there should be a "bond" 
formation between the two layers. In our case, we want 
to show that this bond formation is basically due to the 
presence of /-electrons of Dy through a direct (or indi- 
rect) interaction. Indeed, we have seen that if we con- 
sider the / electrons in the core (i.e. using the VASP 
code, we use the Dy_3 potential), we don't observe any 
exchange striction effect. If we consider them in the va- 
lence (and so let them eventually interact), we do have 
an exchange striction effect, i.e. the two ferromagnetic 
sheets approach each other. Although the effect is small, 
it disappears completely when using Dy_3. This is clearly 
confirmed by the following computational experiment: 
starting from the FEi ionic structure, we freeze the / 
electrons in the core (using Dy_3 POTCAR file within 
VASP). Obviously the Dy atoms are not spin-polarized in 
this case. We let the system to relax to its new electronic 
and ionic ground state. Not surprisingly, we found that 
the system relax to a non-polar state (PE state). Vicev- 
ersa, starting from the PE ionic structure, and treating 
the / electrons in the valence, the FE state is stabilized. 
In summary, when considering the / electrons as valence 
states, the PE state becomes unstable, the point 
group symmetry is spontaneously broken to C211 and the 
system evolves towards a stable and polar state. If the / 
electrons are removed from the valence and frozen in the 
core, the PE state remains stable. This unambiguously 
confirms that / states are a necessary ingredient for fer- 
roelcctricity in DyFeOs. So, coming back to our ELF 
function, we consider the difference in ELF (DELE) be- 
tween the situation when / electrons are in the valence 
(which stabilizes the FE state, as previously observed) 
and when they are frozen in the core (thus stabilizing 
the PE state), for the same ionic configuration (for in- 
stance the CS), i.e. DELF(r)=ELF/„„, (f)-ELF/^„^^(r). 
The physical interpretation is as follows: positive values 
of DELE show up in regions where the electron local- 
ization is higher, i.e. the bonding between FM layers 
is strengthened. It can be showniSl that a positive iso- 
surface of DELE projected into the ab plane is mainly 
localized between FM layers and, more specifically, in 
the region between Dy and O],'^ . This points therefore 
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to a bonding interaction between FM layers mediated by 
Ot,t_ 

^ eq 

Our electronic structure analysis interpret this find- 
ing as an efficient mediation of 0-2p and Dy-c? between 
the relevant Dy-4/ and Fe-Sd electrons. For details see 
Ref J^. These results pave therefore the way to the in- 
teraction between / and d electronic states as an addi- 
tional degree of freedom to tailor ferroelectric and mag- 
netic properties in multiferroic compounds. Addition- 
ally, the fact that the Dy and Fe interaction is mediated 
by 0-2p states suggests possible routes to tailor the FE 
polarization. For instance, compressive or tensile strain 
along the polar axis might change the octahedral tilting, 
favoring or disfavoring the interaction via the interme- 
diate O states, eventually leading to a change in the FE 
polarization. From the methodological point of view, Dy- 
ferrite constituted a nice benchmark for the theoretical 
treatment of 4/ electrons, usually a hard task within ab- 
initio approaches. Our results were shown to be robust 
with respect to the different state-of-the-art computa- 
tional schemes used for d and / localized states, such as 
the DFT-I-U method, the Heyd-Scuseria-Ernherof (HSE) 
hybrid functional, and the GW approachiSi. 



2. tt4-i spin arrangement in CdV2 04 

Spinels exhibit several unusual features, which make 
them an emerging class of materials in several fields, 
including magnetoelectricity and multiferroicity. Most 
oxide-based spinels contain trivalent and divalent cations 
and have the general formula A+^B+'^204. The structure 
type is that of the MgAl2 04 mineral, which is cubic and 
contains eight formula units. Oxygen ions show a cubic 
close-packed arrangement and the cations occupy both 
octahedral and tetrahedral sites. When the tetrahedral 
sites are occupied by divalent cations only and the octa- 
hedral sites by trivalent ions only, the structure is termed 
a normal spinel. 

The series of Vanadium oxide spinels with A = Cd, 
Mn, Zn or Mg and B=V is particularly interesting, since 
it approaches a Mott transition when the V-V distance 
is reduced sufficiently, either by applying pressure or by 
changing the size of the A cationii^i^^. ZnV204 is the 
member of the series which is closest to the metallic 
state. The ground state of the spinel compounds has 
stimulated an intense theoretical research in the last few 
yearsii^^-ii^. A tetragonal distortion induces the t2g lev- 
els to split into a lower dxy level and a twofold degenerate 
dxz and dyz level. It is clear that the first electron of 
(d^) occupies the dxy level, whereas the second one is lo- 
cated in a combination of the other t2g-orbitals (dxz and 

Three main models have been proposed to describe 
the dxz and dyz occupations: the "real" orbital-ordered 
model, where the ground state consists of alternating oc- 
cupation of the dyz or dxz orbital in adjacent layers along 
the c axis; the "complex" orbital order model, where the 



second electron occupies the complex orbital (dyzidxz) 
which has an unquenched value of the orbital angular mo- 
mentum; a third model takes into account the proximity 
of ZnV2 04 to the itinerant-electron boundary, showing a 
dimerization along the V-V chains, which is described 
by the formation of homopolar molecular V-V bonds, 
characterized by a partial electronic delocalization (the 
orbital wave functions in this case would be a real com- 
bination of orbitals (dyz^dxz) with no net orbital mo- 
ment). As far as the magnetic structure is concerned, 
chains in the xy (ab) plane show an AFM ordering, while 
chains in the xz and yz plane have t-lt-l' spins. The over- 
all spin structure is thus collinear. It must be noted that 
magnetic A sites like Mn or Fe have a much more compli- 
cated magnetic structure. The multiferroic behaviour in 
the spinel class of materials is very rare. The only few ex- 
ceptions are CoCr204, HgCr2S4 or CdCr2S4. Usually, a 
more complex magnetic structure is involved there, such 
as spiral magnetism. 

A multiferroic behaviour was recently reported in a 
ternary spinel, CdV2 04 (hereafter called CVO) with 
collinear antiferromagnetic ground state, where the ferro- 
electricity arises from a local-exchange striction mecha- 
nism. This shows that not only spiral magnetism can give 
rise to polarization in spinels, but also collinear structures 
through exchange striction, thus broadening the class of 
multiferroic systems by including this class of oxides. 




FIG. 4. Perspective view of the spin and crystal structure of 
the CdV2 04. 

In Figure S] we show a perspective view of the unit 
cell of CVO, along with the spin structure. It is easy 
to recognize the ttii chains. In this particular case, 
ab-initio calculations proved to be very useful to high- 
light the microscopic mechanism leading to a finite P 
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especially through a series of computational experiments 
and a trend study as a function of the U parameter. 
Our study can be summarized as follows: 

• We started from the centric lAi/amd space group 
and we imposed tttt spins. In this initial ionic con- 
figuration, there was no V-V dimerization. Then, 
for any U between and 8 eV, we let the ionic de- 
grees of freedoms relax. What about the final ionic 
and electronic ground states? We found: i) no V- 
V dimerization; ii) no inversion symmetry-breaking 
and thus, no polarization (for U sufficiently large 
to keep the system insulating, i.e. larger than U —A 
eV); 

• Starting from the centric lAi/amd symmetry, we 
then imposed ttii spins along the [101] and [Oil] 
directions. After ionic relaxations, we found: i) 
the formation of short (S) and long (L) bonds be- 
tween tt and tJ- spins, respectively; ii) inversion 
symmetry breaking and appearance of a finite po- 
larization. Although we found it difficult to unam- 
biguously extract the final symmetry group of the 
relaxed structure, due to numerical noise, a com- 
patible symmetry group may be the space group 
80, C4-6, Ml (The threshold on the atomic position 
for the symmetry check has been fixed to 0.001 A) 



a) Centric phase 



b) Polar phase 




U-J £bV) 



FIG. 5. Sketch of the V spin chain in the (a) centric and 
in the (b) polar phase. Lower panel: Relevant V-O-V angles 
(as defined in b)) along the the spin chain as a function of 
the effective Coulomb parameter (U-J) within the DFT-I-U 
formafism. (Adapted from Fig. 1 of Ref.-^^.) 

Further inspection in the mechanism for ferroelectric- 
ity, shows that: the ttii spin order, imposed onto the 



centric lAi/amd structure, gives rise to an electronic in- 
stability that ultimately results in i) a V-V dimerization 
and ii) formation of short and long V-0 bonds, compat- 
ible with a staggered xz,yz orbital ordering. It is worth- 
while to note that this electronic driving force towards i) 
and ii) shows up already before performing ionic relax- 
ations: the two t (4-) V sites are inequivalent and such 
inequi valency, upon relaxations, drives the V-V dimer- 
ization; the two oxygens bonded to t,t V (or 1,1 V) are 
inequivalent, in turn giving rise, upon ionic relaxations, 
to a weakly staggered orbital ordering. We believe that 
both effects, i.e. dimerization and orbital ordering, co- 
operate to induce polarization. In Figure [5] we show the 
centric phase with a FM spin chain, and the polar phase 
with an tt4-4- spin chain. In the centric phase, all V-O-V 
angles are equivalent along the chain, see Figure [S] (a) . 
Note, however, that V-0 distances are slightly different, 
due to the peculiar coordination of the spinel structure: 
Each O is an "apical" one with respect to one V ion and 
a "planar" one with respect to the other neighboring V 
ion. This is in principle compatible with the presence of 
partial orbital ordering, even in the FM spin chain. As 
expected from the centrosymmetric space group, no po- 
larization is found from our calculations for this case. In 
the polar phase, see Figure[S](b), (i) the angles ai and /?2 
(q!2 and /?! ) become inequivalent due to the formation of 
short and long V-V bonds; (ii) ai and /3i become differ- 
ent. The long V-0 bonds are compatible with a weakly 
staggered xz,yz orbital ordering. As a result, local dipole 
moments, originating from the inequivalency of oxygens, 
appear due to different ai and /32 {0-2 and /3i) angles; 
since the dipoles do not compensate, we observe a net P 
in the unit cell. Further details can be found in RefJ^. 



C. Charge ordering 

Transition metal oxides often show correlated elec- 
trons, which, under certain conditions (due to a complex 
interplay among Coulomb repulsion, electron-phonon in- 
teraction, Jahn- Teller effects, etc), can lead to electronic 
charges being localized on different ionic sites (the lat- 
ter belonging to the same chemical species) in an or- 
dered fashion: this phenomenon is labeled as charge- 
disproportionation or Charge ordering (CO) and it is a 
(first- or second-order) phase transition, with well de- 
fined critical temperatures. Charge ordering (CO) was 
proposed as a phenomenon that can induce ferroelectric- 
ity in those cases where, similar to spin ordering, the 
symmetry of the CO pattern below the critical order- 
ing temperature lacks inversion symmetryi^^. Several 
examples were put forward: Polar CO was identified in 
Fe^+/Fe3+ occurring in magnetito^'^'^^^i^^", Mn3+/Mn''+ 
in half-doped manganite o^^^'^'^^ , Ni^+/Ni''+ in rare-earth 
nickelates— . Ab-initio calculations showed that a po- 
lar CO can lead to potentially large polarization (of the 
order of few /iC/cm^), so that it can be considered as 
an efficient mechanism in the context of electronic ferro- 
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electricity. Among the mentioned systems, however, one 
has to make an important distinction as for the driving 
mechanisms: there are cases, such as nickelates and man- 
ganites, where multiferroicity shows up when both spin 
and charge orders occur. It is actually their combination 
that drives ferroelectricity (see Figure [6] a). In this situ- 
ation, magnetism has a relevant role in the development 
of the dipolar order (actually, magnetism and ferroelec- 
tricity share the same physical origin) and a large ME 
coupling is therefore expected. In other cases, such as 
magnetite, deeply discussed in the following paragraph, 
magnetism seems not to be involved in the ferroelectric 
transition. In this large magnetoelectric coupling 

is a-priori not to be expected, also given the different 
origin and ordering temperatures for the magnetic and 
dipolar orders. 



1. Fe^'^ /Fe^'^ charge patterns in magnetite 
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FIG. 6. a) Coexistence of charge and spin orders, whose co- 
operation breaks inversion symmetry and gives rise to fer- 
roelectricity. Shown is the case of RNiOs (R — rare earth), 
where Ni - nominally trivalent - charge-disproportionates into 
Ni^"*" and Ni^^ and magnetically orders inan^ — t^J-^i 
fashion along the [111] pseudo-cubic direction. Due to the Ni 
spin and charge inequivalence, a dimerization is induced along 
[111] and polarization arises, b) Pure charge-ordering giving 
rise to polarization. Shown is the case of Fe in magnetite 
along the b direction, where a dimerized chain of Fe^"*'-Fe^^ 
forms, due to a complex interplay between Coulomb repul- 
sion in the Fe tetrahedral network, entropy, electron-phonon 
interaction, etc. Black thin arrows denote Fe displacements 
with respect to the centrosymmetric configuration. 



FIG. 7. Ionic structure of Fe octahedral sites in a) P2/c and 
b) Cc. Orange and blue balls show Fe^^ and Fe'^^ ions, re- 
spectively. Fe tetrahedra of 2:2 and 3:1 CO patterns are high- 
lighted by yellow and black color planes, respectively. Electric 
dipole moments caused by charge shifts are indicated by red 
arrows. For details, see Ref.-^^. In the four squared blue 
boxes, we show the charge/orbital ordering of Fe minority t2g 
states in PE and FE states at different planes for different z 
internal coordinates in the unit cell: (upper-left, al): PE for 
z=3/8; (lower-left, a2): PE for z=2/8; (upper-right, bl): FE 
for z=3/8; (lower- right, b2): FE for z=3/8; . 

Magnetite is probably one of the (if not "the") most 
studied magnets: it was discovered in Greece around 
6*'* century before Christ and, since then, it has al- 
ways attracted lots of interests. Fe304 (formally Fe^'*' 
jPg2.5-(-pg2.5+j^ Q2--J g]^Q.^^g g^j^ inverted cubic spinel 
structure with a Fd — 3m space group at room temper- 
ature. The inverted spinel structure shows Fe^ and Fes 
ion sites coordinated to O ions, i.e. tetrahedral Fe sites 
are occupied by Fe^ ions, whereas octahedral Fe sites 
are occupied by Fes ions. The latter form a network 
of corner sharing tetrahedra. Magnetic moments on Fe^ 
sites are antiparallel to those of Fes sites, so that ferri- 
magnetism is the ground state. Magnetite undergoes a 
first order metal-insulator transition (called Verwey tran- 
sition) at around 120 Ki2^, where the resistivity increases 
by two orders of magnitude. Correspondingly, the crystal 
structure changes from cubic to monoclinic. Verwey has 
proposed the metal-insulator transition to originate from 
charge ordering at Fe^ sites (Fe^+ [Fe^+Fe^+]B Ol^). 
The pattern of the charge ordering, however, constitutes 
a matter of debate and it is still unknown. Anderson 
has pointed ou1>i2^ that, when putting two Fe^"*" and two 
Fe'^"'" sites on each Fes tetrahedron (so called "2:2" pat- 
tern), the number of Fe^^-Fe^^ ion pairs is maximized , 
therefore giving rise to the lowest possible energy from 
the Coulomb repulsion point of view. However, the An- 
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derson criterion for low temperature magnetite is incon- 
sistent with recent experimental results and alternative 
patterns with "3:1" CO arrangement (three Fe^"*" and one 
Fe'^"'" ions in a tetrahedron, or viceversa) or even "mixed 
75% 3:1 and 25% 2:2" have been put forward. 

In our recent works, pure charge-order (CO) was care- 
fully investigated as a potential source of inversion- sym- 
metry-breaking electronic order. Indeed, this mecha- 
nism was explored in magnetite below the Verwey metal- 
insulator transition. In a joint theory-experiment studj*^^ 
and in following purely theoretical studie9i22,ii2£, we 
showed that magnetite in the Cc symmetry (predicted 
by density-functional-theory (DFT) to be the ground 
state and suggested to experimentally occur at very 
low temperatures) shows a non-centrosymmetric CO of 
Fe^^/Fe'^"'" on octahedral Fes sites of Fe304, with P ~ 
5 /iC/cm^. Magnetite might therefore be considered as 
"one of the first multiferroics known to mankind" . Re- 
markably, it is a beautiful example from another point of 
view: since what is usually searched for in electrically- 
controllable spintronic devices is a net magnetization, 
magnetite, being a Jerri-magnet, overcomes the limita- 
tion of a zero (and therefore uncontrollable) magnetiza- 
tion that occurs in many other multifcrroic antiferromag- 
nets. 

As shown in Figure [71 octahedral Fe sites, arranged 
in the corner sharing tetrahedral network are located in 
xy planes with z=i/S (i—Q...l). The P2/c paraelectric 
state has C26 -f (0, 0, 1/2), /, aab + (0, 0, 1/2) symme- 
tries (with a full 3:1 tetrahedron CO arrangement) and 
the Cc ferroelectric state has E, (J2b + (0, 0, 1/2) symme- 
tries (with a mixed CO pattern) in a conventional base- 
centered monoclinic cell so that there are two equivalent 
atoms (cfr B12 and B12' sites in Figure [7]). We remark 
that the mirror symmetry along with the translation vec- 
tor forbids any net polarization along b and finite P is 
allowed only along the a and c directions. The difference 
between the two Cc ferroelectric and P2/c paraelectric 
CO distributions (see Figure [7] a) and b)) can be under- 
stood when assuming a charge "shift" from B12 to B14 
site and in the upper part of the cell, from B12' to B14', 
all the other sites keeping their valence state unaltered. 
Each "charge shift" creates two 2:2 CO tetrahedra, so 
as to form, in total, four 2:2 tetrahedra in the unit cell 
(cfr. Figure [7]). The resulting CO pattern lacks inver- 
sion symmetry, therefore allowing FE polarization. The 
Berry phase approach predicts quite a large polarization, 
its direction lying in the ac-mirror plane. The DFT re- 
sults are in excellent agreement with recently reported 
experimental values for magnetite thin films (reporting 
P of the order of 5.5 /iC/cm^ in the ah plane with the c 
component not measured) as well as with earlier experi- 
ments on single crystals {Pa — 4.8 /iC/cm^ and Pc — 1.5 
fxC/cw?). We've verified that the polarization values are 
not largely affected by the value of the Hubbard U pa- 
rameter, as shown in Table ID What we also note is that, 
upon increasing U and keeping the atomic configuration 
fixed to that obtained for J7=4.5 eV, the charge separa- 



tion between Fe^"*" and Fe^"*" is increased, in agreement 
with what intuitively expected: a "full charge dispro- 
portionation" to occur in the limit of an infinitely large 
Coulomb repulsion. 



TABLE I. Charge separation (cs, i.e. difference of d-charges 
between Fe'^''" and Fe'^'*' ions in the atomic sphere with lA 
radius) and the corresponding FE V Berry (in /iC/cm^) vs 
Coulomb repulsion U (J is fixed to 0.89 eV). 

U (eV) 4.5 6i) 8.0 

cs 0.17 0.23 0.30 

^Berry (-4.41, 0, 4.12) (-4.42, 0, 4.81) (-4.33, 0, 5.07) 



We remark that, especially in CO-materials (such as 
magnetite), we' ve often found an antiferroelectric (AFE) 
phase energetically competing with a ferroelectric on o^^° . 
For example, in magnetite, the ground-state FE Cc sym- 
metry is only a few meV/unit-cell lower than AFE P2/c; 
also in the case of Fe304 phases in which an interme- 
diate bond-and-site-centered CO occurs, the AFE P2/c 
symmetry competes with the FE P2 symmetry 



D. Orbital ordering 

In addition to charge and spin, electronic degrees of 
freedom include the orbital one. Indeed, many transition 
metal oxides clearly show, below a critical temperature, 
an orbital-order (00), often driven by the Jahn- Teller 
effect and accompanied by structural distortions in the 
octahedral or tetrahedral oxygen cages which surround 
transition metal ions. In principle, there seems to be 
no obstacle, from the symmetry point of view, to the 
fact that 00 itself could break inversion symmetry and 
give rise to polarization, i.e. nothing precludes orbital- 
induced ferroelectricity. However, there are no estab- 
lished examples where this happens in a clear and simple 
way. For example, the Ruddlesden-Popper bilayer man- 
ganite, 

Pr(Sro.i Cao.9)2Mn207J^, was proposed as a candidate 
material in this context, since the rotation of orbital pat- 
tern below a defined critical temperature was found to 
happen along with a ferroelectric state. However, many 
degrees of freedom were active at the same time in that 
compound: in addition to 00 and ferroelectricity, CO 
coupled with the underlying lattice distortion was also 
occurring, so that the link between 00 and polarization 
is actually under debate. A pure system in which 00 
by itself drives ferroelectricity is still to be found. What 
we will discuss in the following section is, on the other 
hand, a compound where orbital-order occurs and, indi- 
rectly via hydrogen bonds, induces a polar state in an 
organic-inorganic hybrid^. 
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1. Cu-based metal- organic- framework: role of orbital order 
and hydrogen-bond 

Metal-organic frameworks (MOFs) are an interesting 
new class of materials made up of extended ordered 
networks of metal cations linked by organic bridgeai^. 
They are hybrid organic-inorganic materials at the in- 
terface between molecular chemistry and materials sci- 
ence. There is an huge interest in these materials for their 
potential technological applications such as gas storage, 
exchange or separation, catalysis, drug delivery, optics, 
magnetismi^ii^. Furthermore, due to their dual na- 
ture, they can be engineered in almost infinite ways by 
playing with the organic/inorganic componentsi^^ii^. 
A very recent family of MOFs, with a more dense 
topology, mimics the ABX3 perovskite inorganic topol- 
ogy. These compounds show interesting magnetic, op- 
tical, electronic and dielectric properties. Last but not 
least, coexistence of ferroelectricity and magnetism, i.e. 
multiferroicityi^^-i^. 




FIG. 8. Crystal structure of the Cu-MOF: (a) side view, (b) 
a single linear chain along the c polar axis, (c) an octahedron 
with the HCOO organic linkers, (d) perspective view of the 
Guanidinium ion. 

We have recently studied a Cu based MOF, namely 
[C(NH2)3]Cu[(HCOO)3], first synthesized by Ke-Li Hu 
et alJ^ In this compound, hereafter called Cu-MOF, A 
is the guanidinium ion C(NH2)"'", B is the Jahn- Teller 
Cu+'^ ion with d^ with tjC^ electronic configuration and 
X is the carboxylic linker HCCO". At low temperature, 
it crystallizes in a polar space group Pna2i. Further- 
more, a magnetic study revealed that it displays spin- 
canted antiferromagnetism, with a Neel temperature of 
4.6 K. In addition to the general spin-canted antiferro- 
magnetism, it is a magnetic system with low dimensional 
character>i^. A magnetic ordering in a polar space group 
immediately calls for a possible niultiferroic behaviour, 
although no ferroelectric hysteresis loop has been mea- 
sured yet. Our theoretical study predicts and supports a 



multiferroic behaviour. Furthermore it highlights inter- 
esting features in this appealing class of materials, such 
as an unusual microscopic mechanism for ferroelectricity 
and the magnetoelectric effect. 

In Figure [SJ we show: (a) side view of the Cu-MOF; 
(b) a single chain of octahedra connected by the HCOO 
organic linkers along the polar c axis; (c) Cu octahedra 
with HCOO groups; (d) perspective view of Guanidinium 
ion. The Jahn- Teller distortion of the Cu octahedra give 
rise to two short (s) and two long (1) equatorial Cu-Oe<j 
bonds with lengths ^ 2.0 A and 2.4 A respectively, and 
two medium (m) apical Cu-Oap bonds. The cooperative 
Jahn- Teller distortion is characterized by CuOg octahe- 
dra elongated along the [1,1,0] and [1,1,0] in the ab plane. 

Our calculations show that the most stable mag- 
netic configuration is the AFM-A type, which shows 
ab intra-plane ferromagnetically aligned spins which are 
inter-plane antiferromagnetically coupled. When in- 
cluding spin-orbit coupling in the calculation, a weak- 
ferromagnetic (FM) component arise due to a small cant- 
ing of the spins. The weak-FM component Ma is along 
the a axis, perpendicular to c axis. The presence of 
weak-ferromagnetism is in agreement with experimental 
observationi^. The perovskite Cu-MOF is very simi- 
lar to KCuFs which is considered as a prototypical sys- 
tem for a cooperative Jahn- Teller effect, orbital ordering, 
and a quasi-one-dimensional antiferromagnetic Heisen- 
berg chain. In fact, the Cu-MOF shows a particular type 
of orbital order, in which a single hole alternately oc- 
cupies 3da;2_^2 and 3dj^2_^2 orbital states of the Cu+^ 
ions (3d^ electronic configuration)^^. The cooperative 
Jahn- Teller distortion is characterized by Cu(HC00)6 
octahedra alternatively elongated along the perpendicu- 
lar [1,1,0] or [1,1,0] directions in the a&-plane, i.e. giv- 
ing rise to an antiferrodistortive pattern. It is important 
to note that the anti-ferro-distortive (AFD) modes are 
usually non-polar distortions in standard inorganic per- 
ovskite like compounds, and, as such, they should not 
give rise to ferroelectric polarization. Despite this, our 
calculations show that the AFD distortions in Cu-MOF 
are strictly correlated to the presence of the polarization. 

Our study highlights very interesting properties of this 
compound: 

it is ferroelectric with an estimated polarization P of 
0.37 fxC/crri^, with polar axis along c; 
ii) the microscopic mechanism is very intriguing: we 
found that non-polar AFD distortions are intimately re- 
lated to ferroelectric polarization suggesting that they 
may be the "source" of ferroelectricity; 
Hi) inspection into the microscopic mechanism of ii) 
shows that AFD distortions acting on the BX3 frame- 
work are coupled to A-group ions through inter- 
vening hydrogen-bonds between the Oxygens of the 
Cu(HC00)6 and the H atoms of the A-group. While 
the AFD distortions alone would preserve centrosymme- 
try, the O • • • H bonds induce asymmetric distortions into 
the A-group, ultimately responsible of the presence of 
dipoles mainly localized at the A-group which, in turn. 
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give rise to a finite polarization; 

iv) the weak- ferromagnetic component is strictly corre- 
lated with the ferroelectric polarization: when P is equal 
to zero, the weak-FM components goes to zero; at the +P 
state, it is +Ma and at the ~P state it is —Ma- There- 
fore, our calculations predicts that the Cu-MOF should 
be a magnetoelectric multiferroic: it should be possible 
to control the magnetization by an external electric field. 
In particular, an electric field along the c axis, which 
would switch the spontaneous polarization, would at the 
same time, switch the sign of Ma ■ Although both P and 
M are small in magnitude, this opens new avenues in the 
multiferroic research in such a novel and exciting class 
of materials. It goes without saying that there could 
be large room for engineering these compounds for en- 
hancing these magnetoelectric effects due to the organic- 
inorganic duality characteristic of MOFs. More discus- 
sions about the Cu-MOF can be found in Ref 

In conclusion, MOFs are materials at the border- 
line between chemistry and solid state physics and MF- 
MOFs represents a "dual-bridge" between the two fields, 
exploiting knowledges from the inorganic as well as 
organic material science. We expect that this dual- 
bridge will be the source of new and interesting phys- 
ical propertiesi^ii^, which ab-initio studies can easily 
unveil^^. Incidentally, we note that ab-initio character- 
ization of MF-MOFs are almost totally lacking in the 
current literature, and our recent study^ is certainly en- 
couraging in terms of interesting results. 



E. Charge-spin dimers in donor-acceptor TTF-CA 

In comparison with inorganic materials, organic com- 
pounds have been synthesized in large number but fer- 
roelectric properties have been found only rarely in that 
class of materialsi^. A breakthrough in organic ferro- 
electricity was recently achieved by the discovery of very 
large room-temperature ferroelectric polarization in the 
croconic acid, a well-known low-molecular- weight organic 
compound22. It is believed that it may be fruitful to 
search among known - but poorly characterized - organic 
compounds for organic ferroelectrics with enhanced polar 
properties suitable for device applications^^ . 

Coexistence of ferroelectric and magnetic order in or- 
ganic materials is an even rarer property. Recently, we 
have predicted by ab-initio calculations that multifer- 
roicity may be found in TTF-CA molecular crystal^^, 
a multi-component molecular system which produces a 
typical displacive-type ferroelectricity by displacing op- 
positely charged species. TTF-CA is a charge-transfer 
(CT) complex composed of electron donor (D) and accep- 
tor (A) molecules, such as tetrathiafulvalene (TTF) and 
tetrachloro-pbenzoquinone (CA)i^".i^. This compound 
is particularly interesting because it shows a neutral-ionic 
(NI) phase transition, i. e. a transition between a van der 
Waals molecular assembly to an ionic solidi^ii^. The 
ionized molecules form DA dimers, D+^'-A"'', where p is 



a degree of charge transfer, with a lowering of the crystal 
structure from Pn to a polar P2i/n space group, where 
the originally non polar D---A---D---A sequence with 
regular intermolecular separation are symmetry broken 
to a polar chain formed by the DA dimers characterized 
by the formation of pairs of short and long bonds along 
the stacking axis a. In essence, above the NI transition 
temperature Tni of ^ 84 K^^^ the system is in a neutral 
and paraelectric state with D+''- • • A^''- ■ • D+''- ■ • A^'' 
with p=0.2-0.3. Below T^j the system becomes fer- 
roelectric with a stacking D+''- • • A"'' D+''- ■ • A"'' 

with p'^ 0.6, with an ionic and ferroelectric state charac- 
terized by a Pierls-like dimerization. 

In this framework, we have performed ab-initio cal- 
culations by using the recently introduced screened hy- 
brid functional Heyd-Scuseria-Ernzerhof (HSE)^^. The 
use of the hybrid functional has been particularly impor- 
tant here for several reasons. First, it is important to 
improve the description of the HOMO-LUMO gap which 
governs the degree of charge transfer between molecu- 
lar units; second, we found it impossible to stabilize a 
magnetic state by using the local or semilocal approxi- 
mation to the exchange-correlation functional (LDA or 
GGA); third, the commonly used DFT+U method for 
improving the electronic structure of "strongly" corre- 
lated system can not be directly applied here: the reason 
is that in molecules, the localized orbitals are multicenter 
rather than single-center, since the basis set correspond 
to ortho-normal molecular orbitals instead of orthonor- 
mal atomic orbitals. The main results of our study can 
be summarized as follows: the HSE ground state of the 
TTF-CA crystal shows an antiferromagnetic (AFM) or- 
dering, more stable than a non-magnetic one by ~ 80 
meV per unit cell. Note that starting with an initial fer- 
romagnetic (FM) configuration, the solution converges 
again to an AFM one. This demonstrates the robust- 
ness of our AFM solution. The single TTF and CA 
units become spin-polarized with a net polarization of 
^ 0.40 pb per molecule. The HSE energy gap of spin- 
polarized dimerized state is 0.5 eV. A ferroelectric state 
with a coexisting magnetic ordering as ground state char- 
acterizes the TTF-CA as the first multiferroic organic 
crystals, predicted by ab-initio calculations. The calcu- 
lated ferroelectric polarization in the AFM state is 3.5 
pC/cw? while in the NM state is 8.0 pC/cm^. The 
sensitivity of the polarization to the magnetic state is 
mainly due due to the large increase of the electronic 
component of polarization upon changing from the NM 
to AFM state. The ionic one, on the other hand, does 
not depend much on the magnetic state, and it is al- 
ways opposite to the electronic one. Further details can 
be found in the original article^-. Finally, the multifer- 
roicity in TTF-CA has been independently confirmed by 
theoretical calculations^^. A recent experimental result 
show that one-dimensional quantum magnets, such as or- 
ganic charge-transfer complexes, could be promising can- 
didates in the development of magnetically controllable 
ferroelectric materialsi^"-i2£. 
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F. Coupled distortions in NaLaMnWOe 

As well discussed so far, perovskite (with formula 
ABX3) is one of the crystalline structures which is most 
commonly occurring and most important in all of mate- 
rials science. Because of the great flexibility inherent in 
the perovskite structure, mainly due to the corner shar- 
ing octahedra, there are many different types of distor- 
tions which occur starting from the ideal cubic structure. 
These include tilting of the octahedra, displacements of 
the cations out of the centers of their coordination poly- 
hedra and distortions of the octahedra driven by elec- 
tronic factors {i.e. Jahn- Teller distortions). Many of the 
physical properties of perovskites depend crucially on the 
details of these distortions, particularly the electronic, 
magnetic and dielectric properties which are so impor- 
tant for many of the applications of perovskite materials. 

The new class of double perovskites AA'BB'Og intro- 
duces yet another degree of freedom, namely the possibil- 
ity of cation ordering on both A and B sites. Obviously, 
this greatly increases the possibility of functional design 
in this class of compounds. More than 20 new examples 
of this structure type have been discovered so far. These 
materials are found to have highly complex microstruc- 
tures and show potential for multiferroic behaviori2ir-i2&. 

We have recently presented a theoretical study of the 
structural and ferroelectric properties of the new double- 
perovkite NaLaMnWOgi^^-iS, which belongs to this 
class of materials, by combining group-theoretical anal- 
ysis and first-principles calculations to explore the ori- 
gin of the polar state in this compound. NaLaMnWOe 
orders magnetically at low temperature in a polar space 
group. However, the ferroelectricity has neither been cal- 
culated nor measured yet. We found that ferroelectricity 
originates not from a usual type of lattice distortion in- 
volving small off-centerings of ions, as usually occurring 
in the prototypical ferroelectric BaTiOs, but from the 
combination of two oxygen rotational distortions. 

This idea of rotation driven ferroelectricity is a very 
exciting recent development in the field of ferroelectrics 
as well as in the related field of multiferroics. This rep- 
resents an interesting new route to produce new mul- 
tiferroic and magnetoelectric materials, relying on the 
idea of starting with non-polar materials and then in- 
duce multiple non-polar instabilities; under appropriate 
circumstances, this can induce a ferroelectric polariza- 
tion, as first predicted in Ref based on general group 
theory arguments and analyzed in the SrBi2Nb209 com- 
pound by means of a symmetry analysis combined with 
density-functional theory calculations by Perez-Mato et 
alJ^. In that case, ferroelectricity was found to arise 
from the interplay of several degrees of freedom, not all of 
them associated with unstable or nearly-unstable modes. 
In particular, a coupling between polarization and two 
octahedral-rotation modes was invoked to explain the 
behavior^i^. Bousquet et al. have demonstrated that 
ferroelectricity is produced by local rotational modes in 
a SrTiOa/PbTiOa superlattice^. Benedek and Fennie 



proposed that the combination of two lattice rotations, 
neither of which produces ferroelectric properties individ- 
ually, can induce a ME coupling, weak ferromagnetism, 
and ferroelectricity^^^. Indeed, we now know that ro- 
tations of the oxygen octahedra, in combination^— "-i^^ 
and even individuall y^^^'^^^ , can produce ferroelectricity, 
modify the magnetic order, and favor magnetoelectricity. 

Our study on NaLaMnWOe represents another step 
forward along this new emerging direction. In Figure 
ini (a) we show a perspective (a) and side (b) view of 
the magnetic unit cell of the compound. In refJ^ we 




FIG. 9. Crystal structure of the NaLaMnWOs compound: 
(a) perspective view and (b) side view. The spins are shown 
as arrows only in (a). 

showed that this compound is a potentially very inter- 
esting multiferroic compounds for several reasons: i) the 
estimated polarization is very large, about 16 /zC/cm^; 
ii) an intriguing mechanism is at the basis of the fer- 
roelecticity: two primary non-polar distortions such as 
tilting and rotation of octahedra - typical of perovskite 
systems - in combination with cation ordering induce the 
breaking of inversion symmetry and allows for a ferro- 
electric polarization. By comparing the low symmetry 
structure with a parent phase of P4/nmm symmetry, two 
distortion modes are found dominant. They correspond 
to MnOe and WOg octahedron tilt modes, often found 
in many simple perovskites. While in the latter these 
common tilting instabilities yield non-polar phases, in 
NaLaMnWOg the additional presence of the A-A cation 
ordering is sufficient to make these rigid unit modes as 
a source of improper ferroelectricity. Through a trilinear 
coupling with the two unstable tilting modes, a polar dis- 
tortion is induced: a negligible polar instability does ex- 
ist, but the additional A cation layer ordering makes fer- 
roelectrically active some tilting modes of the octahedra 
that in simple perovskites and in B-ordered double per- 
ovskites only give rise to non-polar behaviour. Despite 
its secondary character, this polarization is coupled with 
the dominant tilting modes and its switching is bound to 
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produce the switching of one of two tilts, enhancing in 
this way a possible interaction with the magnetic order- 
ing. Through a trilinear coupling with the two unstable 
tilting modes, a significant polarization is induced. We 
hope that this study will stimulate further investi- gation 
of cation ordering as a tool to convert ubiquitous well- 
known steric non-polar instabilities into mechanisms for 
producing improper ferroelectrics, as well as new multi- 
ferroics. Further details can be found in Ref 



IV. CONCLUSIONS AND PERSPECTIVES 

It is evident that multiferroics currently represent a 
vivid and promising field of materials science, the en- 
thusiasm being mainly driven by i) the impression that 
multiferroics are much more common than what origi- 
nally thought and ii) the rich and still largely unexplored 
variety of mechanisms leading to the coexistence of mul- 
tiple orders. Density functional theory is able to identify 
the microscopic mechanisms, their strengths as well as 
their limitations, their chemical and physical origin, so it 
seems a particularly suited technique for the analysis of 
this complex class of materials. 

As shown in this paper, our recent activity was devoted 
to those materials where ferroelectricity is induced by 
peculiar charge, spin or orbital orders which lack inver- 
sion symmetry, i.e. to the so called "improper electronic 
ferroelectrics" . We'll try in this conclusive paragraph to 
summarize the main findings of our recent activity as well 
as to give guidelines towards an efficient materials-design 
for optimized multiferroics. 

• When dealing with spin-driven ferroelectricity, two 
mechanisms have been mainly explored so far: the 
first one is based on Heisenberg exchange coupling, 
whereas the second one is based on relativistic 
Dzyaloshinskii-Moriya interaction. An important 
issue regards the efficiency of these mechanisms, or 
equivalently, the magnitude of the polarization in- 
duced by the two mechanisms. According to our es- 
timates based on density functional theory, the po- 
larization caused by Heisenberg exchange (at play 
in collinear spin configurations) is much larger than 
that driven by Dzyaloshinskii-Moriya exchange (at 
play in non-collinear spin configurations) . This was 
quantitatively shown, for example, in nickelates^, 
when we compared first -principles estimates ob- 
tained for polarization in both collinear or spiral 
spin arrangements, as experimentally proposed^^, 
the two values differing by approximately two or- 
ders of magnitude. This is consistent with what 
expected, based on the argument that relativistic 
effects are not very sizeable in 3ci transition metal 
oxides and that the symmetric exchange coupling is 
much larger than the corresponding antisymmetric 
component. 

• Many systems were studied, all of them showing 



Heisenberg-driven polarizationiiiSliS^. The mag- 
nitude of the latter, as expected, is strongly de- 
pendent on the involved transition metal and spin- 
state. In a rather qualitative, general and naive 
way, one can invoke the size of spin moments and 
exchange-interactions to rationalize the different 
behaviour in different oxides. In particular, we have 
quantitatively estimated the biggest effect to occur 
in rare-earth manganites involving Mn^"*" (d^), with 
polarization of the order of jjiC /crn?. This big value 
is reasonably due to the large size of the spin mo- 
ment ('^ 4 ^b) and to rather strong interactions 
between Mn Cg and oxygen p states which are able 
to well mediate the Mn-Mn exchange interaction in 
ortho-manganites. As a result, the (large) generic 
Heisenberg term, J^^ Si ■ Sj, can induce apprecia- 
ble changes in the energy and related sizeable (lo- 
cal) distortions, depending on whether Si and Sj 
spins are parallel or antiparallel. A large polar- 
ization is therefore expected when all the (local) 
distortions are summed up over the magnetically- 
ordered unit cell, given the overall polar distor- 
tion pattern and non-centrosymmetric spin config- 
uration. For a similar reason (but producing an 
opposite result), the V-V dimerization is able to 
induce a much smaller polarization (order of few 
tenths of a jiC/cm^) in spinel Cd-based vanadate 
(cfr Sec lIII B"2]) : spin moments are smaller in size 
and the states at play are t2g {i.e. much less prone 
than Eg to interactions with oxygen, due to their 
main non-bonding character, with likely smaller ex- 
change constants and related smaller distortions) . 
Also in the case of spin-driven ferroelectricity for 
f — d systems, of which the prototypical DyFeOa 
case was discussed in Section lHI B 11 one can expect 
a smaller exchange interaction between localized 4/ 
states and semi-localized^ 3d states, compared to 
the exchange interaction between 3c? states. Since 
it is the f — d coupling that induces the polar con- 
figuration in DyFeOa and the spins have to order 
on both Dy and Fe sublattices, we expect a smaller 
polarization than in, say, orthomanganites, as in- 
deed predicted by our ab-initio simulations. Inci- 
dentally, we note that spin-driven ferroelectricity 
based on 4/ states is likely to develop only at tem- 
peratures as low as < 10 K, where the rare-earth 
ions order, as is the case of DyFeOa 

• How to increase the ordering temperatures for 
magnetically-induced ferroelectrics (commonly of 
the order of few tens of Kelvins) represents in gen- 
eral one of the toughest challenges towards find- 
ing a so-called "killer-app" in the field of mul- 
tiferroics for them to become really technologi- 
cally appealing. Although progresses were made 
in recent years (for example, by focusing on com- 
bined charge- and spin-ordered materials, such as 
nickelates*^ and hole-doped manganites show- 
ing ordering temperatures of 150-200 K), room 
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temperature operation is still a dream. One of the 
limitations might be constituted by the fact that 
electronic magnetic ferroelectrics are generally frus- 
trated materials (showing either spin or charge or 
orbital frustration), whose ordering temperatures 
are - due to competing interactions - intrinsically 
small. The possible way out could be to deal with 
large exchange-coupling constants, which is what 
happens in nickelates and in CuO tenoriteiSSiiSi. 
The latter material, in particular, offers an exam- 
ple of spin-spiral-based multiferroic with a large or- 
dering temperature (230 K). 

• When dealing with CO-induced ferroelectric- 
ity, polar CO (as identified in Fe^+/Fe^+ oc- 
curring in magnetitc^^^, Mn'^4-/Mn^+ in half- 
doped manganites^'^^, Ni^+/Ni''+ in rare-earth 
nickelates^) can lead to potentially large ferroelec- 
tric polarization (of the order of few nC/cm?). In 
this respect, CO-driven ferroelectricity is an effi- 
cient mechanism and more work towards a bet- 
ter understanding should be definitely carried out. 
There is an important distinction to be made as 
for the microscopic origin of ferroelectricity. In 
nickelates and manganites, multiferroicity shows up 
when both spin and charge orders occur: it is actu- 
ally their combination that drives ferroelectricity. 
In this case, magnetism has a relevant role in the 
development of the dipolar order (actually, mag- 
netism and ferroelectricity share the same physical 
origin) and a large magnetoelectric (ME) coupling 
is therefore expected. In magnetite, on the other 
hand, where CO alone drives ferroelectricity, mag- 
netism seems not to be involved in the ferroelectric 
transition (Fe304 becomes ferrimagnetic at around 
860 K and remains so down to very low temper- 
atures, all over the metal-insulator Verwey transi- 
tion) . In this large magnetoelectric coupling 
is a-priori not to be expected, also given the differ- 
ent origin and ordering temperatures for the mag- 
netic and dipolar orders. A possible reason for the 
magnetoelectric coupling to occur (as reported in 
a very early study back into 1994^^^, involves the 
relativistic spin-orbit coupling in the peculiar Cc 
symmetryiS^ 

• As a general rule-of-thumb for a larger CO- 
induced polarization, we remark that, in the 
case of non-centrosymmetric CO, the charge dis- 
proportionation (CD) should be maximized. In 
fact, according to a picture based on point-charge 
dipoles, valid mostly for systems where the bond 
is largely ionic, this would guarantee a larger 
polarization. To our experience, we focused on 
oxides (magnetite, Fe304) ^^i-^^^i^^° and fluorides 
(Ko.6Feg+Fei^+F3)iM, both with Ye^+ /Ye^+ charge 
disproportionation: the comparison between mag- 
netite and fluorides showed that a larger CD oc- 
curs when iron is bonded to Fluorine rather than 



to Oxygen. Therefore, choosing a more ionic com- 
pound seems promising to achieve a large CD and 
related higher polarization. 

• On the theory-side, as from the methodological 
point of view, we remark that the treatment of cor- 
relation effects is often important to get a quan- 
titatively reliable description of multiferroics. In 
addition to the common DFT+U approach, that 
we have used in a variety of studies^i^^ii^, we 
have carried out careful simulations of prototypi- 
cal multiferroics, such as the "proper" BiFeOa and 
the "improper" AFM-E HoMnOs, using a state-of- 
the-art hybrid exchange-correlation functional, ob- 
tained by mixing the non-local Fock-exchange with 
a "standard" parametrized exchange-functional, 
with encouraging results^. In particular, 
hybrid-functionals with the "ideal" mixing between 
3/4 local and 1/4 non-local exchange-correlaton 
potential, appear to give an accurate description of 
structural, electronic, ferroelectric, magnetic prop- 
erties for most of the studied materials. This 
same technique was later applied to the TTF- 
CA organic^^ and DyFeOa multiferroicaiSl. A 
correct description of the vibrational as well as 
spin-phonon coupling within hybrid functionals was 
also recently shown for many well-characterized 
oxideai^l. 

• We have abundantly shown that ferroelectricity 
driven by electronic degrees of freedom can occur 
(and has actually been experimentally observed) 
in many systems where spin-order and charge- 
orders drive the rising of polarization. On the 
other hand, ferroelectricty induced by orbital or- 
der has remained for long elusive. In Sec lIII D"T] 
and in Ref.— , we focused on a class of materi- 
als called metal-organic frameworks, i.e. corre- 
sponding organic-inorganic hybrids of perovskite 
crystals. This architecture is much more flexi- 
ble (due to organic groups instead of single oxy- 
gen anions) and chemically more rich than usual 
inorganic perovskites (in terms, for example, of 
organic polar or non-polar groups which occupy 
the empty site corresponding to the A-site cation) . 
This chemical richness suggests that new mecha- 
nisms might arise in this class of materials. In 
particular, we considered a MOF based on Cu^''" 
ions at the center of octahedral cages of COOH- 
groups and with guanidium molecules occupying 
A sites. A delicate interplay between Jahn- Teller 
distortions around Cu^+ (in turn related to the 
antiferrodistortive orbital-ordering) and hydrogen 
bonding with guanidinium groups induces a small 
ferroelectric polarization. Moreover, following a tri- 
linear coupling between magnetization, antiferro- 
magnetism and polarization allowed by symmetry 
in the ferroelectric crystal, we reported a linear pro- 
portionality between weak-magnetization (induced 
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by Dzyaloshinskii-Moriya coupling) and polariza- 
tion, pointing to the long-sought electrical control 
of magnetization. 

• Most of the niultiferroics discovered so far are an- 
tiferromagnets (due to usual strong superexchange 
in oxides which often favors antiparallel spins), so 
that their technological appeal is poor. In order 
to overcome this limitation, we mention two pos- 
sible solutions: i) choose a compound where fer- 
rimagnetism (quite a common spin configuration, 
occurring in spinels, such as magnetite, or in dou- 
ble perovskites, etc) is the magnetic ground state, 
so as to show a net magnetization that can be well 
controlled via a magnetic field, ii) consider (multi- 
ferroic)-antiferromagnets exchange-linked to ferro- 
magnets, so as to build an artificial heterostructure 
where both electric and magnetic degrees of free- 
dom are simultaneously activ&i2i"-i2S. There, the 
phenomenon of exchange-bias can be used, for ex- 
ample, to control the magnetization of a FM over- 
layer by means of an electric field which primarily 
modifies the ferroelectric as well as the antiferro- 
magnetic properties of a multiferroic layer, which 
the FM overlayer is, in turn, exchange-coupled to 
(proposals for applications in this direction already 
came®). 

• The TTF-CA donor-acceptor organic crystal was 
probably one of the first examples of organic fer- 
roelectrics treated from first-principles^® (inciden- 
tally, we remark that a breakthrough in the field 
was later achieved in 2010^, when the croconic acid 
in crystalline form was discovered to be ferroelectric 
with large polarization persisting at least up to 400 
K, showing an excellent qualitative and quantita- 
tive agreement between ab-initio theory and exper- 
iments). When looking at TTF-CA, the combina- 
tion of charge transfer and structural dimerization 
results in an opposite behavior for the electronic 
and ionic contribution to polarization, which was 
first predicted from first-principles and later con- 
firmed by experiments. The field of organic crystals 
might be richer of ferroelectrics than what origi- 
nally thought, and efforts should be devoted in the 
near future to investigating polarization in many 
charge- transfer salts, charge-ordered systems and 
other strongly-correlated organics. Proposals to- 
wards this direction already appeared in the litera- 
ture, for example pointing to the quasi-two dimen- 
sional organic sah a-(BEDT-TTF)2l3i^. 

• Among the many different routes to new mul- 
tiferroics, particular interest has been raised by 
the idea of a trilinear coupling among polariza- 
tion and different octahedral distortions. The 
idea has been originally proposed for ferroelec- 
tric Aurivillius compounds^®^ and recently redis- 
covered in the context of layered manganites, such 



as Ca3Mn207>i2^. In the present work, we've sug- 
gested yet another possibility of trilinear coupling 
in double-perovskites with formula AA'BB'Og, 
based again on functional octahedral distortions 
as well as cation ordering and resulting in a large 
polarizatioE^iS^. 

The field has to face several challenges in the coming 
years. While electronic ferroelectrics show at least two 
characteristics which are definitely appealing for techno- 
logical applications (i.e. expected giant magnetoelectric 
coupling and expected ultrafast switching), the real bot- 
tleneck is represented by the fact that working temper- 
atures of the most studied "improper" multiferroics are 
too small for applications and efforts should be devoted 
to the increase of operating temperature range. 

An alternative way - much closer to applications - 
might be represented by interfaces between prototypical 
ferromagnets and prototypical ferroelectrics - both with 
high operating temperatures - combined in artificial het- 
erostructures where the magnetoelectric coupling could 
be engineered and optimized. This vivid field has not 
been treated in the present review, which was mostly de- 
voted to "bulk" multiferroics; however, it should be kept 
in mind that many progresses were made in recent years 
on this kind of artificial systems (such as Fe/BaTiOa^SS), 
aimed in the end at implementing multiferroic memories. 

Possibilities to overcome the present limitations in bulk 
electronic ferroelectrics might involve either the discov- 
ery of new physical mechanisms in known materials or 
the optimization of known mechanisms (such as spin or 
charge-order induced ferroelectricity found for transition- 
metal oxides) in "new" materials (such as organic or 
organic- inorganic hybrids, novel complex oxides). In 
general, a better understanding of spin-phonon cou- 
pling in transition metal oxides through first-principles 
approaches might certainly help in designing materials 
with large polarization and strong magnetism (possibly 
ferromagnetism)^^"*. In this respect, some works have re- 
cently investigated the possibility of ferroelectric and fer- 
romagnetic instabilities in transition-metal-oxides with 
cations under strain or volume expansion (i.e. Ca- 
based manganitesi^iiSSir^SS^, or La-based cromites22i) . 

In summary, the physics of electronic ferroelectrics 
is rich and complex, so that surprises in both mecha- 
nisms and materials are to be expected in the coming 
years . In this respect, we remark that efficient and re- 
liable modeling approaches can greatly contribute to the 
field, by proposing new materials, new mechanisms and 
their quantitative estimates. We therefore hope that the 
present manuscript will contribute to stimulate further 
scientific interests from the experimental point of view 
towards this peculiar class of materials . 
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